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The slow-rotation approximation of Hartle is developed to a setting where a charged rotating 
fluid is present. The linearized Einstein-Maxwell equations are solved on the background of the 
Reissner-Nordstrom space-time in the exterior electrovacuum region. The theory is put to action 
for the charged generalization of the Wahlquist solution found by Garcia. The Garcia solution 
is transformed to coordinates suitable for the matching and expanded in powers of the angular 
velocity. The two domains are then matched along the zero pressure surface using the Darmois- 
Israel procedure. We prove a theorem to the effect that the exterior region is asymptotically flat if 
and only if the parameter C2, characterizing the magnitude of an external magnetic field, vanishes. 
We obtain the form of the constant C2 for the Garcia solution. We conjecture that the Garcia metric 
cannot be matched to an asymptotically flat exterior electrovacuum region even to hrst order in the 
angular velocity. This conjecture is supported by a high precision numerical analysis. 
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I. INTRODUCTION 

There are a surprizingly few rotating perfect fluid so- 
lutions of Einstein's relativistic equations known to date. 
One of the most comprehensive of these is the rigidly ro- 
tating charged perfect fluid solution of Garcia Q]. This 
metric is type D in the Petrov classification of the curva- 
ture. The Einstein-Maxwell equations are satisfied with 
a stress tensor which is the sum of that of the perfect 
fluid and the Maxwell field. The fluid medium carries 
electric charge, and it has a divergence-free four-current. 

Because of the tantalizing lack of explicit, relativisti- 
cally rotating fluid stellar models in the literature, it is 
only natural to ask if it is possible to join a suitably cho- 
sen domain of the Garcia solution to an external source- 
free Einstein-Maxwell region. 

The Garcia solution reduces to the Wahlquist space- 
time in the absence of electric charge. Unfortunately, 
as we have proven in an earlier publication the 
Wahlquist solution is unsuitable as a model of an isolated 
relativistic object. A priori it is unclear if the presence 
of electric charge can supply the necessary ingredients for 
a smooth matching. 

The purpose of the present paper is to search for the 
conditions of matching the Garcia metric to an asymp- 
totically well-behaved external electrovacuum domain. A 
suitable vehicle for performing this research is the slow- 
rotation approximation. We could successfully use this 
approach in our previous computations for the Wahlquist 
space-time The foundations of the slow- rotation ap- 
proximation were laid down by Hartle |^] . In its original 
form, the Hartle method matches an internal fluid do- 
main with an external vacuum. An essential tenet of this 
approach is the independence of the first-order perturba- 
tion function 10 of the polar angle For a generic setting 



in the presence of the Maxwell field, one would have to go 
through the proof of this independence. In the present 
case, however, we can omit this general theory since 10 
is a function of the radius alone in the slowly rotating 
Garcia space-time. As a consequence of the matching 
conditions, then, this property is inherited by the elec- 
trovacuum exterior. 

In this paper, we apply the slow-rotation approxima- 
tion to the Garcia space-time. Our result is that the elec- 
tric charge works against the conditions of matching. We 
find that matching to an asymptotically flat electrovac- 
uum exterior is impossible already to the first order in 
the angular velocity. This comes as a surprize to us since 
all uncharged fluids can be matched to an asymptotically 
flat vacuum exterior to first order in the rotation param- 
eter. This is because such vacuum space-times do not 
differ from the first-order form of the the Kerr metric. 

The paper is organized as follows. In the next section, 
we compute the Garcia metric in various forms that will 
be necessary to launch the matching process. In sub- 
section 



[I B we take the static limit which, of course, is 

we as- 



spherically symmetric. Next, in subsection II C 



scmble the slowly rotating form of the metric. In Sec. 
Ill we prepare the field quantities in the external elec- 



trovacuum domain. Thus we get the general solution for 
the first-order perturbation function uj, and we investi- 
gate its asymptotic properties. The actual matching of 



the two domains takes place in Sec. IV. As before, we re- 
quire that the Darmois-Israel conditions H are satisfied. 



The details are worked out in subsection IVA] , general- 
izing the theory in the presence of a Maxwell field. The 
solution of the matching conditions for the static state 
and to first order in t he an gular velocity is carried out in 
subsections IV B and IV C , respectively. The results are 
further discussed in Sec. |V|. The Appendix is devoted to 
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the issue of the choice of independent parameters. 



are 



II. THE INTERIOR METRIC 

The charged generalization of the Wahlquist solution 
given by Garcia M has the metric 



ds 2 



^ (dr + SNdaf - ^ (dr + SMdaf 



A 



d£_ 
~P 



df_ 

Q 



(i) 



with the real functions 



A 

AI 

N 



M-N 



1 



sinh 2 (fca;) — £q 



1 



(2) 



sin 2 (fcy)-Co 
P = a + h [ 2n + x ( a + P 2 )] sinh(2fcr) 
+ b - (g + fix) 2 cosh(2fcir) 

Q = a ~ Jk [ 2m + v( a + P 2 )] sin(2fcy) 
b+ (e + f3y) 2 ] cos(2fcy) . 



Here a, b, e, g, k, m, n, a, /?, 8 and £o are eleven con- 
stants. The number of independent parameters modulo 
diffeomorphisms is eight ( Cf. Appendix) . This metric is a 
solution of Einstein's equations G a fj — 8n + T^J 

with 

where p is the pressure, /j, the density, u a the velocity of 
the fluid, and F a p is the electromagnetic field tensor. In 
comoving coordinates, the velocity vector is u — u d/dr. 
The pressure and density are given by 



8*p = -— (Q- P) + ak 2 + T, 
k 2 

8tt^ = 3— (Q - P) - ak 2 - £ 



E = 



2(3k 



[(e + 0y) sm(2fcy) + (g + (3x) sinh(2fca;)] 



The electromagnetic four-potential A a determines the 
field tensorf] as F a p — d a Ap — dpA a , and its components 



*Garcfa Q uses the definition F a p — 2A[ at p], which gives 
rise to an overall sign difference in A a . 



A X Ay — 



A T 



A a = -7-rr [(e + (3y) Msin(2ky) + (g + fix) N sinh(2kx)} 



2Ak 



[(e + fiy) sin(2fcy) + (g + f3x) sinh(2fca;)] 



4k 2 /3 ' 

The electric current vector j a — pu a satisfies the field 
equation 47rj Q = The current Airj a = 2k 2 f35? has 

no divergence (as it should, since the Maxwell equations 
are satisfied) and the conserved charge density is 



k 2 P Q-P 



2tt V A 

Note that this solution has been named the Wahlquist- 
Newman metric in p]. 

A. The Wahlquist form 

Our goal in this section is to bring the metric (Q) to 
a form which reduces to that of the Wahlquist solution 
in the no-charge limit. The static limit of the Wahlquist 
form will, in turn, yield the charged generalization of the 
Whittaker space-time. To achieve this, we transform to 
the new coordinates £ and £ by setting 



fc£ = smh(kx) , k( = sm(ky). 
Then we have 

a = e 2 + c 2 ■ 



(3) 



(4) 



In order to be able to introduce Wahlquist's functions hi 
and h 2l we rescale the functions P and Q using a new 
constant parameter ro; 

hi = r 2 Q , h 2 = r 2 P. 



To choose the factors of the other terms appropriately, 
we also rescale the r and a coordinates as 



1 



1 

t=—T. 



The metric takes the form 



hi 



hih 2 



ds 2 = — _ Ad(f f + r 2 (5 2 A 

til — Il2 



A 



dip 2 



Ar 2 



d( 2 



(i-k 2 e)hi (i + k 2 e)h 2 

where A is Wahlquist's function, 

'hie+h 2 e 



Sr 



hi - h 2 



(5) 



(G) 



2 



We introduce the rescaled constants 



m = mr n 



e = er 



(7) 



9 = gr , (3 = (3r Q 



and three new constants as follows. In place of a we shall 
use the constant re, the constant a will be replaced by C 
and finally E will take the place of the constant b, using 
the following definitions^ 



- 2 =ri (a + f3 2 
C = t$(a + b) 



(8) 



E 



-2r 2 6fc 2 - - 2 . 



We get for the functions 



, c 

hi = —y 



C v/l — fc 2 C 2 arcsin(fc£) 

k 



ro 







e + — arcsin(fcO 
fc 



1 - 2k 2 C) +C + EC 



(9) 



/l 2 = -4r 



£ v/l + fc 2 £ 2 arcsinh(k£) 

k 



ro 



<7 + ^arcsinh(k£) 
fc 



(1 + 2k 2 i 2 ) +C-E? 



(10) 



The electromagnetic potential 1-form A = A a dx a is 



.4 



1 

'A 



arcsinh(k£) 



x [dt + (C 2 + Co) Sr dkp] 



1 

A 



_ /3 

e + — arcsin(fcC) 
fc 



x [dt - (£ 2 - £o) *r dv] • 
The pressure and density become 



(11) 



87175 



(hi - h 2 ) + 4-2 (l-/? 2 « 2 )+S (12) 



Ar 2 



fc 2 



A: 2 



8?r/i = 3 {hl _/,,)_ (!_ £V) _ E (13) 

^ r r K 



(14) 



The expressions (p|)-(p^|) contain a large number of pa- 
rameters: C, -B, m, 6, fc, re, e, g, ,9, ro, 5 and £o- How many 



^We have slight differences from the notation of Garcia Q: 
m and b have different factors, k includes a term /3 2 , and _E 
and C are not set to 1 yet 



of these twelve constants are necessary to uniquely de- 
scribe the metric? This question will next be addressed. 

By using linear coordinate transformations of t and (p 
we can obviously set 6 = 1 and £0 = 0. The parame- 
ter ro has been introduced in order to enable ourselves 
to go to the slow-rotation limit. The scaling of this pa- 
rameter can be chosen arbitrarily. Defining r' = cro, 



C" 



2 C, E' = c 2 E, fh 



1 - ^3, 



y - ~3 



c 3 b, k' = k/c, 



e' = ce, g' = eg, /?' = c/3, where c is a constant, the dC, 2 
and d£ 2 terms in the metric remain unchanged with the 
primed quantities. However, we need to perform a trans- 
formation of the coordinates t = ct', <p = c 2 (p' to keep 
the structure of the other components. This shows that 
the two sets of constants {C, E, fh,b,k, re, e, g, j3, ro } and 
{C, E',fh' ,b' , fc, re' , e' ,g' , 0' ,r } provide parametriza- 
tions of a single physical state equivalent under the dif- 
feomorphism {£', tp' , t'} = {(, £, c~ 2 (p, c -1 i}. 

Another freedom of choice follows from the coordinate 
transformation £ = /£, £ = /£, where / is a constant. 
In the metric feh we have 



e 2 + c 2 



1 

T 2 



i' 2 + c" 2 



f 2 



rA 



(15) 



Introducing the constants fc = fc//, C = f 4 C, E = 
f 2 E, m" = prfi, b" = f 3 b, re" = K,/f,e" = f 2 e, g" = 
f 2 g, (3 = f/3, and defining h 1 and h 2 similarly to h\ 
and /12, we get 



In 



1 



/ 



a hi 
4 J- 



1 



4 h 2 



(16) 



Then the dC, 2 and d£ 2 terms transform appropriately. 
However, to get the correct transformation for the other 
terms in the metric, we need to change the coordi- 
nates as t = ft and (p = f 3 Lp . It follows that the 
two sets of constants, {C, E,rh,b, fc, re, e, /3, ro} and 

,m ,b ,fc , re ,e ,g ,/? , r 1 are equivalent de- 
scriptors. 

At this stage, the coordinates are completely fixed. 
Hence there are eight independent physical parameters. 
The corresponding analysis of the parameters in the origi- 
nal coordinates is not needed in the present paper, but for 
convenience, is included as an Appendix. Using the coor- 
dinate transformations x 1 ' = x 11 (x fe ) and x l " = x l " (x k ), 
for example, the constants C and E can always be made 
±1 or 0. Garcia 0] discusses the case C = E = 1, 
which for zero e, g and (3 yields the original uncharged 
Wahlquist solution. When C or E can only be rescaled 
to —1 or 0, some different solutions might arise, as was 
first noted by Mars and Senovilla H in the uncharged 
case. 

There exist two combinations of these transformations 
which preserve many of the constants. Table 1. shows the 
sets of constants giving metrics equivalent to the original 
set {C,E, m, b, k, re, e, g, 0, ro} , when choosing cf = 1 
and cf 2 = 1 . 



3 



and 



TABLE I. Equivalent sets of parameters 











1 


C E fh 


b k 


k e g (3 r 


c/=l 


c 2 C E fh 


b i* 

c 


K ce eg (3 ±r 


c/ 2 = l 


C f 2 E ffh 


pb fk 


jk e g f/3 f 2 r Q 



B. The static charged fluid 

In this section we consider the static limiting case 
of the Garcia space-time. Following the procedure of 
Wahlquist fjj, we introduce a new radial coordinate 

r = (r (17) 

and a new constant 7 in place of k by substituting 

k = 77-0 (18) 

everywhere. The constant 7 is related to Wahlquist's p a 
by 7 2 = n 2 p s . The static limit can be obtained then by 
going to zero with rg. In the limit tq — > 0, 



lim -±hi 

r ->Q \ r 2 - 



where 



hi=E- 
1 

+ — 



2m 



\f\ — 7 2 r 2 



(19) 



(20) 



1 \f\ — 7 2 r 2 arcsin(7r) 

7r 

2 



arcsin(7r) 



(1 - 2 7 V) 



r 

and e = roe. The function /12 has the limiting form 

r-2 



h 2 = lim hi 



C-££ 2 -26£-(g + ^) (21) 



where b = b/r^. Using the limits r 2 A — ■> r 2 , and A — > 0, 
the metric becomes 



ds z 



-hidt 2 



dr 2 



(1 - 7 2 r 2 )/ ll 



r 2 [^+5 2 ~h 2 d^ 



(22) 



A simpler form arises if we introduce a new radial co- 
ordinate z by setting 



Then 



"fr = sin z. (23) 



hi = E — 2m7 cot 2 H — =■ (1 — z cot z) 

AC 

+ (ij + fiz) 2 (cot 2 z- 1) (24) 



ds 2 = -/ix^ 2 



dz 2 sin z ( d£ 



S 2 h 2 dip 2 



(25) 



The only nonvanishing component of the electromagnetic 
potential is 

A t = — (£7 + /3z) cot 2 . 
The pressure @ and the density @ arc 



87rp = 7 2 I -/ll + - /?' 

8tt^ = 7 2 f 3Ai - 7^ + /3 2 
E = 4/3 7 2 A 4 . 



-E 
E 



(26) 



The bracketed term in (25) is the metric of a two 



surface. We next show that there is a parameter range 
for which this two-surface is the two-sphere S 2 . In order 
to get the metric of the two-sphere, we introduce a new 
coordinate 9 in place of £ by S 2 h 2 = c\ sin 2 9, where ci is 
some constant. Then 



2 dh 2 d9 
— - = 2ci sin 9 cos 9— 



di 



(27) 



and 



d£_ 2 _ 4c 2 sin 2 (9 cos 2 6> / dh 2 
17 ~ S*Ka V^T 



d9 2 



4 c x - S 2 h : 



The condition 



S 2 



de 



dh 2 



d9 2 



(28) 



= c 2 d9 2 



(29) 



where c 2 is another constant, is satisfied if and only if 
1 



c 2 = 



E + [3 2 



(30) 



d = S 2 



C-g 2 + 



E + P 2 



We can use the map in the third row of Table 1 to set 
c 2 = 1 by making E + /3 2 = 1. Of course, this can be done 
only when E + p 2 is positive. The Lorentzian signature 
of the metric (123) hold whenever it is possible to set 



(b + gp)' 

E + P 2 



(31) 



4 



using the second row of Table 1 . Finally we set S = 1 by 
rescaling the tp coordinate. With these choices h,2 = 1 — 

(l- + b + g/3j and we end up with the simple coordinate 
transformation 



Then 



£ + b + g/3 = cos ( 



hi = 1 — (3 2 — 2to7 cot z H — - (1 — z cot z) 



(32) 



+ {erf + (3z) 2 (cot 2 z - l) 
and the metric becomes 



(33) 



ds 2 = -/iicft 2 + 4t 

7 2 



dz 2 

17 



- sin 2 z i 



sin 2 edip 2 



(34) 



The center, determined by z = 0, can be regular only 
if fh — and e = 0. 



C. The combined transformation 

In the previous sections, we have been led to trans- 
forming the Garcia metric to the Wahlquist coordinates 
and hence to a form which is amenable for accessing the 
static limit. We now combine these procedures into a 
direct transformation, sidestepping the Wahlquist form, 
and compute the quantities for slow rotation. 

The coordinate transformations (||), (17) and (23), 
with k — "fro, can be combined to the single transfor- 
mation 



z = ky . 



(35) 



Garcfa's parameters are expressed in terms of ours as 
follows, 



k = 


jr 


TO 

m = 


n = 


b 
'() 




e = 


e 
r o 


9 

9 = — 

ro 


P = 


I 


C 

a = + 
r o 


b = 


1 


My 


2r^7 2 



(36) 



E 



We now substitute the new parameters given by ( p6[ ) into 
the metric form (Q) and introduce new coordinates by 

y = — , t = r t , a = -r\<p . (37) 
7 r o 



Using Table 1, we find that the coordinate system is fixed 
by setting 



E = 1 - P 2 

2 



(38) 



8=1 , Co = 0, 
C = 1 + f - (b + . 

The regularity at the center is ensured by 

to = , e = 0. (39) 
The angular coordinate $ is introduced by writing 

x = cos-d - b - g/3 . (40) 



The coordinate $ only equals the 9 defined in Sec. [IB 
when the fluid is static. Then for small ro to linear order 
we obtain the metric 



ds 2 = -hidt 2 



dz 2 
l 2 ~hi 



sin 2 z 



d$ 2 + sin 2 (dip - codty 



(41) 



with 

hi = 1 
and 



/3 2 + — (1 - z cot z) + f3 2 z 2 (cot 2 z - l) (42) 



w = r Q 



■ 2 

sin z 



(hi -I). 



(43) 



The first-order calculation shows that uj does not depend 
on the angular coordinate. To this order, the only com- 
ponent of the four-potential A to pick up a small new 
contribution is A v : 

A = —j3z cot z dt 

+ r a [/3sin 2 i?(l - zcot z) - (3 - 3 cos??] dip. (44) 

Note that to first order in ro, the metric is independent 
of the magnetic monopole charge parameter g. In fact, 
the monopole contribution affects only the Maxwell equa- 
tions but does not affect the gravitational equations. 



III. THE ELECTROVACUUM EXTERIOR 

The metric of the ambient electrovacuum domain, to 
first order in the angular velocity, has the Reissner- 
Nordstrom form modified by the non-diagonal rotation 
term, 



ds 2 



„ to e \ , 9 
1 - 2 — + — ] dt 2 

r r A 



1-2 



+ r 2 dir +sin*i9 (dtp - udt)" 



dr 
(45) 



5 



The four-potential A — A a dx a , where A a = A a (r, 
has a time component 



A, 



(46) 



From the (r,r) and (r, i?) components of Einstein's Eqs. 
we get no contribution to At, and the (t, i?) component 
gives A rj ,j = A^. r . Hence there is a gauge in which 



A. = A* = 0. 



(47) 



The other perturbation components are obtained in the 
remainder of this section. Before doing that, we remark 
that for the special case of the slowly rotating Kerr- 
Newman solution the forms (Q) and ( |47| ) of the four- 
potential components remain valid and we have 



A(KN) 



,(KN) _ 



ae sin $ 



2am ae 2 



(48) 



This is a separable equation. Introducing the separation 
constant K, we get 



sin 2 tf 



(53) 



with the solution 



liCcos 2 ^ Q 



4 



C 5 



cosfl 



(54) 



2 e e c 

We next consider the radial part of Eq. (|52|) . 
(i) For the case e 2 to 2 the general solution of the 
radial equation is 

e 2 - 2mr ^ 2 _ 
^=Ci+ „ s , C + -C 2 
3m J r« r 



e 2 (r 2 — 2 mr + c 2 ) [r + to + (e 2 — r 2 ) I/] 



m 2 (e 2 — m 2 ) r 4 



where 



, , 2e 2 



Co = & 2 + 3eWC 2 + 



C 3 (55) 



(56) 



We now proceed to computing the first-order contri- 
butions for a general electrovacuum. From (ip, t) compo- 
nent of Einstein's equations the second-order differential 
equation follows: 



dr 2 



o dui 
+ 4r 3 — 
dr 



4e 31 



r 2 $ <9r 



0. 



(49) 



For the uncharged case, e = 0, this has the solution u> — 
ujQ + 2am/r 3 where a is the Kerr parameter. The value of 
the constant ojq can be set to zero by the transformation 
ip — » <p + C\t. 

The solution of (^) for a generic charge e ^ is 



A v = ^r 4 sin 2 ^ + /(tf). 
4e dr 

The Maxwell equation for the component A v is 



(50) 



(2 



2 2\ A<n n . 2 n 

m r — r — e ) + r sin t/ e — — 

or ar 



(51) 



+ 2 I m 

r 



9r 



^ ^cot^ = 0. 



dti 2 



Substituting A v from (|50|) : 



(2to r 



2\ 4 

e r 



dr 3 



(52) 



+ (| rm „2r 2 -|e 2 ) r 3 — 



(4m-|r) r 3 



dr 



sin 2 



and 
L(r) 



r— yn+-\/rn, : -' — e '-' 



2 \/m 2 — e 3 r — m — \/ m 2 — e 3 



1 / 7T 

v 2 



arctan 



Vc a — m £ 



with the derivative 

dL (r) 



if m > e 
if to < e 
(57) 

(58) 



(ii) For the equilibrium case e 2 = to 2 , the radial solu- 
tion is 



2 5r 2 — 4mr + to 2 



(59) 



15 m(r — m) 2 r A 

When an asymptotically nonrotating frame is chosen, we 
have that Ci = 0, and we shall assume this to hold in 
the sequel. 

Comparing with the Kerr-Newman metric with rota- 
tion parameter a, we obtain 

K = —3am. 

Substituting the solution (j5|), (|^) and in the 
potential (p0|) we get for the case e 2 ^ to 2 : 



1 



ae „ 3e 2 TOr — tot 3 — 2e 4 „ , „ 

— sm 2 H C 2 sm 2 <d 

r 2emr 

3m 2 r(m + r) — 2e 2 m 2 — 4e 4 



6to 3 t (e 2 — to 

3e 2 rnr — tot 3 — 2e 4 

2m 2 r (e 2 — to 2 ) 2 

3am C4 C5 

— 1 1 cos V. 

2e e e 



C 3 esin 2 tf (60) 



C 3 eL (r) sin 2 1? 



6 



For e = ±771 , we get identical limits from both values of 
L in ([57]). The limiting form of the potential A v can be 
obtained by substituting fl59| ) in Eq. (BO): 



am 2 (r + 2m) (r — m) 

: sm v H 

r 2mr 

(2r — rn) (2m 2 — bmr + 5r 2 ) 
15m 2 r (r — m) 3 

3a C4 C5 

— H 1 COSTA 

2mm 



C 2 sin 2 t9 (61) 



C 3 sin 2 7? 



The constant C4 is inessential since it does not appear in 
the Maxwell tensor. 



Theorem. 

When C 2 has a nonzero value, the external solution is 
not asymptotically flat. 

It is therefore important to investigate the values of 
C2 that the matching provides. This investigation will 
be carried out in the next section. 

The term containing C5 is the potential of a magnetic 
monopole. In the given frame, it gives rise to a purely 
radial magnetic field which has the asymptotic form ~ 
1/r 2 . 



IV. MATCHING 



In order to clarify the role of the constants C2 and C5 , 
it will be helpful to introduce a tetrad with the compo- 
nents 



ei 



(-?, 







m 




2 


— + 






r 






m 


-0- 


2 


— + 






r 


1 S 






r di) 






1 




S 



-1/2 



d_ 

dt 



(62) 



1/2 



r sin 77 dtp 



Or 



■sin 1? 1-2 



St' 



In this tetrad the leading terms for large r in the compo- 
nents of the Riemann tensor are given by 



R 



0113 



R 



Co sin 77 



0223 



(63) 



2^0123 — -^0213 



"-710312 — 



Co COS 77 



As seen the terms do not fall off sufficiently fast for the 
space-time to be asymptotically flat. Thus this object 
does not seem to be isolated. To see this we look at the 
asymptotic behaviour of the electromagnetic field in this 
tetrad. From Eq. (pw we obtain 



lim F13 

r — >oo 

lim F 23 



Ci 
e 

C 2 



sinT? 



Bo 



(64) 



COS 77 = B\ . 



This can be interpreted that for nonzero values of C2, 
the fluid ball is immersed in a constant (in this approx- 
imation) external magnetic field, parallel to the axis of 
rotation, and extending to infinity. Note that the Rie- 
mann tensor falls off even though the electromagnetic 
field tends to a constant value. One would expect that 
the contribution of the electromagnetic stresses to the 
curvature will appear in a higher approximation. 

Hence we have the following 



In this section, the matching procedure is described. 
The next subsection outlines the method in general 
terms. We then join the static, spherically symmetric ex- 
ternal and internal domains in subsection IV B The static 



internal state is parameterized by the constants /3, k and 
7. As a result of the matching, these three constants 
determine the radius r\ of the matching surface and the 
parameters m and e of the vacuum exterior. It is a con- 
sequence of the matching conditions that the surface of 
matching coincides with the zero pressure surface of the 
interior. In subsection IV C we carry out the matching 
to first order in the rotation parameter. This will yield 
the parameters of the slowly rotating electrovacuum re- 
gion in terms of rg , the parameter describing the angular 
velocity of the fluid, and in terms of the parameter g. 



A. The matching conditions 

We want a global model from which (i) any surface 
charges and currents, furthermore (ii) surface layers of 
matter are absent. From the first condition it follows 
H that the electromagnetic stress tensor can be contin- 
uously matched at the surface E if we assume that both 
the permeablity and dielectric coefficients are equal in 
the electrovacuum and in the interior. This means with 
(ii) that there is no discontinuity in the pressure p across 
E. 

We write the metric ds 2 = g a pdx a dx^ for both the 
interior and the exterior regions in curvature coordinates, 
{a; Q } = {t, r, 7?, ip} , in the following form: 



ds 2 = ~A 2 dt 2 + B 2 dr 2 



d'd 2 + sin 2 7? {dtp - udt)' 



(65) 



where A, B and to are functions of the radial coordinate 
r alone. Both to order zero and one in the angular ve- 
locity, the constant-pressure surfaces of the perfect fluid 
coincide with the constant r surfaces. 

We match the hypersurface given by r = r\ , of the 
interior region, with the corresponding matching surface 
at r = n , of the exterior region, such that the induced 
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metrics ds 2 \s and induced extrinsic curvatures K\s are 
equal. The continuity of the functions A and to across 
the matching surface can be achieved by transforming 
the coordinates such that 



t = C 6 t', tp = if' + fit' 



(66) 



where Cq and fl are suitably chosen constants and x 1 = r 
and x 2 = •& are unchanged. We shall, however drop the 
primes from the new coordinates. 

The normal of the hypersurface £ has the form 



ld_ 



(67) 



The extrinsic curvature has the nonvanishing compo- 
nents: H 



^00 = 


fflto, in 1 = 


B ' 


(68) 


K 03 = 


5503, l" 1 = 


-^ sin2 ^(^),, 


(69) 


K22 = 


5522a" 1 = 


r 
B 


(70) 


K 33 = 


5S33, in 1 = 


^ sin 2 0. 


(71) 



The junction of if 22 with Eq. ( [70| ) implies that the func- 
tion B must be C at r = r%. (This implies the matching 
of 7^33 .) Next we conclude from Eqs. ( p8| ) and (69), 
respectively, that A and lu are C 1 functions at r\. 

The four-potential A — A a dx a in both regions has a 
small A v component, and A r = A.g = 0. The timelike 
component A t retains its spherically symmetric form to 
first order in the angular velocity. From condition (i), 
the following components of the electromagnetic stress 
tensor F a p are continuous across the matching surface 
£ : F afs h*h p s and F af3 h%<nP ', where h% = 5% - n a n 7 is 
the projection tensor into the orthogonal complement to 
the normal n of the tangent spacefj]. We must match the 
following nonvanishing components: 



F r t 


= At, r 


(72) 


TP 

A rip 


— A 

— n ^,r 


(73) 




= Aiprf. 


(74) 



B. Spherical matching 

We first carry out the matching of the electrovacuum 
region at the sphere r = r% with the internal region at 



II we allow for surface currents and charges (but still assume 
that the permeability and dielectricity are that of the vacuum) 
these conditions take the general form [F a p\ h"hg — and 

[F a p] h^n 13 — 47r jj"^ 1 ' ace ^ where the jump is indicated by a 
bracket. 



z = zx to order zero in the rotation parameter ro- The 
metric of the perfect fluid takes a simpler form when 
using the coordinate z. However, the matching process 
is more transparent when using the radial coordinate r. 
These dual pictures are connected by the transformation 
7r = sinz [Cf. (pSj)]. From the continuity of the metric 
component at the junction surface we find 



n 



sin z\ 

7 



(75) 



Eliminating the mass m between the junction conditions 
of the metric components g tt and K&q, we get the simple 
result for the value of the parameter 



Cq = cos Z\. 



(76) 



Continuity of the radial component of the electromag- 
netic field yields the value of the total charge, 



z\ — sin z\ cos Z\ ) 



(77) 



We next eliminate m between the junction condition of 
gtt and that of the extrinsic curvature component K u - 
Solving for k 2 , we get 



1 tan Z\ 

K 2 Zl 



+ 2f3 2 (2 + zi cot2zi). 



(78) 



The mass is then obtained from the condition of conti- 
nuity of the gtt component as 



m 



n 
T 



COS 2 Zi 



-\ (z? + zi sin cos 2zi + \ sin 2 2zi) . 

27sinzi v i> 

(79) 

The specific charge of the body is characterized by the 
function m 2 — e 2 , given by Eqs. (|77|), ( |78| ) and 
Taylor expanding about the origin, z\ = 0, we get 

A(3 2 



1 



9 7 2 



4 + o (zf) 



(80) 



Hence we find that for small stars, e 2 > m 2 for the values 
|/?| > 1/2. In other cases, we need to treat the different 
types of specific charge separately. 

The pressure at the center is required to be non- 
negative. On Fig. [I], we plot (3 2 as a function of the radius 
zi of the star for the limiting case when the pressure at 
the center vanishes (solid curve). The allowed region of 
p 1 lies under this curve. The other two curves repre- 
sent the two solutions for (3 2 for the extremely charged 
star, m? = e 2 . However, the extremely charged state lies 
outside the physically allowed domain. 
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0.45 




0.25 



FIG. 1. The physically allowed region of the specific charge 
P 2 : The solid curve gives the maximally allowed value as a 
function of the radius. 



C. First-order matching 

From the continuity of the metric component gt v we 
get the angular velocity of the fluid: 



sin z\ 



[/3 2 z 2 (cot 2 zi - 1) - (3 2 + k- 2 (1 - zi cot zx)] 

^( 2mS in Zl - 7 e 2 ) ( 2 m 2 c 2o» n _ 3) 

tri sin z\ cos z\ v A 3( e m J 



_ <= ( e 2 + _ 2m2Hl£i) 



i^i +TO+( e^-2i^£i)i(^) 



47 



sin2zi ^ 2 



(81) 



The condition of continuity of the Maxwell field 
can be investigated by using the form of the four- 
potential in Eqs. ( f44| ) and ([i0|). The matching of the 
magnetic monopolc terms, proportional to cos$, yields 
that 



C 5 



(82) 



The rest of this matching equation, taken together with 
the continuity of F vr and K ttp can be solved for the pa- 
rameters a, C% and C3 as follows, 

ro 



{e/3 _1 7 (4/3 V - 3) - 2/3 2 k 2 zi cos2zi 



(83) 



+2zisin 2 zi +/3k 2 4zi (zi/3 — ej) cot 21 
- (l + 2z 2 ) /? sin2zi} . 



^ ^ „ 33 2 2 -4 -FG-ff r / sin zi 

C2 = Clo - lor ejp z x cos z\ sin zi— — 5— L 

-D 2 V 7 



C 3 



ro /3 4 cos 21 cot zi 
4096 7 4 z 3 e 



FG H 2 . 



(84) 
(85) 



The detailed form of the quantities here is 

/3e7 2 r 



zi tanzi{-512e 4 /3"Vsin b zi 



3L> 2 

- 70 j3 2 zi sin 4 zx sin2;Zi 

+ 4/3 2 z 2 sin 4 zi (39 + 448z 2 + 384z 4 ) 

+ 7/3 2 zi sin 4 zi (5 sin6zi — sin 10 z{) 

+ 2,3 2 z 2 sin 4 zx (6cos2zx + 144 z 2 cos 2z!) 

- 2/3 2 z 2 sin 4 zi (104cos4zi + 9 cos6zi) 
+ 2f3 2 z 2 sinSi (26 cos8zi + 3 cos 10 z x ) 

- 64/3 2 z 3 sin 4 zi (28zi cos4zi + 21sin 3 2zi) 

- 32 /3 2 z 3 sin 4 zi (12 cos 2zi sin 3 2zi + 9 z\ cos %z\ 

- 32/3 2 z 3 sin 4 zi (l28z 2 sin2zi + 24z 2 sin4zi) 

- /3 4 z 2 sin 2 z x (40 - 162 z\ + 1136 z\ + 1536 zf) 
-8^ 4 z 2 sin 2 Z! (z 2 + 56z 4 ) cos2zi 

- /3 4 z 2 sin 2 zx cos 4zi (221 z\ - 1088 z\ - 60) 

- 4 /3 4 z 2 sin 2 zx (6 cos 8 zi — cos 12 z\) 
+ 32/3 4 z 3 sin 2 zx (7 cos2zi - 2) sin 5 2zi 
-I- 4/3 4 z 4 sin 2 zx (3 + 112 z 2 ) cos6zi 

- 3^ 4 z 4 sin 2 zi cos 12 %\ 

+ 256 (3 A z[ sin 2 zi (10 sin2z! +sin4zi) 

+ 128 ^ 4 z x sin 2 zi (5 + 6 cos2zi) sin 3 2zi 

+ 2^ 4 z 4 sin 2 zx [ (31 + 24 z\) cos 8 21 - 2 cos 10 zi] 

+ (3 6 zf B 2 (32 z 3 - 12 zi - 9 sin 2zx + 3 sin 6zx) 

+ 2/3 6 zfB 2 (cos2zi +6 cos4zi -cos6zi) 

+ 128 /^B 2 cos 3 zi sinzi} 



(86) 



where 

D = 16e 2 /3~ 2 7 2 sin 4 z x + /3 4 z\ B 2 

+ 2(3 2 Zl sin 2 z x [ (3 + 24z 2 ) sin 2^ - sin6zi] 
+ 2/3 2 z 2 sin 2 zi (4z x sin 4zi - cos 2zi + 8 cos 4zi ) 
+ 2/3 2 z 2 sin 2 zx (cos 6z x - 8 - 16 zf) 

and 

5 = 4 z 2 + cos 4zi -I- zx sin 4zi - 1 

F = 4/3 2 z 2 + (l + 2/3 2 z 2 ) cos2zi -3/3 2 zi sin2zi - 1 

G = (8 z 2 — l) cos zi + cos 3zi — 4 z 1 sin zx 

H = 2zi (2 + P 2 B) - 4zi cos2zi - 2 sin2zi + sin4zi. 

We next investigate the value of the constant C2 as a 
function of the radius zi and the charge density parame- 
ter (3. The values of C2 are strictly negative in the phys- 
ically allowed region [z x G (0, § ) , (3 G (0, 0.5)} . On Fig. 
^ we display the values of C2 on a 3-dimensional dia- 
gram. For clarity we have chosen the parameter region 
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{zx G (0,0.2) ,p e (0,0.1)} . (The values have been ob- 
tained by the software Maple of the University of Water- 
loo, using a precision of 33 significant digits). 



solve the junction conditions for the electric polarizabil- 
ity P and/or magnetization M . However this would 
probably not give rise to any physically realistic model 
of the substance occupying the interior region. 




FIG. 2. The coefficient C2 as a function of the specific 
charge (3 and the radius Z\. 

Based on this numerical result, we are in position to 
state the following 
Conjecture: 

The Garcia solution cannot be used as the model of an 
isolated rotating body. 

Our conjecture does not prohibit physical applications 
of the Garcia solution. For a nonvanishing parameter C2 , 
the fluid domain is embedded in an external homogeneous 
magnetic field parallel to the axis of rotation. This is a 
typical setting in the interior of galactic disks. 



V. DISCUSSION OF RESULTS 

The Garcia solution is the electrically charged gen- 
eralization of the Wahlquist space time and it carries 
the extra parameter (3 determining the charge density. 
The new degrees of freedom that the presence of elec- 
tric charge bring in the model would seem to raise the 
possibility of a successful matching to an empty exterior 
domain. However, the appearance of these new degrees 
of freedom is compensated by a larger number of match- 
ing conditions. In addition to the surface gravity and 
first curvature, one must match also the electromagnetic 
field, to rule out surface charges and currents. As we 
demonstrated in the main part of this paper, the net ef- 
fect is that the charged generalization of the Wahlquist 
metric is even less likely to serve as a model of an iso- 
lated star. An alternative approach would be to assume 
that the dielectricity or the permeability of the interior 
solution differs from that of the vacuum. We could then 
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VI. APPENDIX 

An important question is how many physical parame- 
ters there are in the metric (|l|). A coordinate transfor- 
mation t = t' + ccr, where c is a constant, can be used to 
set £0 to an arbitrary value. A transformation a = co' 
rescales the constant 5. The determinant of the (r, a) part 
of the metric is ~5 2 PQ, showing that the symmetry axis 
is at those values of x = xq where P(xq) = 0. The coeffi- 
cient of the drda part of the metric is -^(PN — QM). In 
order to obtain a coordinate system which can be made 
regular at the axis we must have M[xq) = 0, i. e. 



£0 



1 



smh(kxo) 



The constant 5 should be set by requiring the usual 
perimeter per radius ratio 2ir for small circles near the 
axis, taking into account the range of the cyclical coor- 
dinate a. 

We should decide whether or not two different sets of 
constants a, b, e, g, fc, m, n, a, /3, <5 and £0 determine 
different spacetimes. Let us assume for the time being 
that we set £0 = and 5=1. If we perform the coordi- 
nate transformation x' = cx, y' = cy and introduce the 
constant k' — k/c then the functions M and N transform 
as 

' ' '^'^ - 1 -sinh^fcV) = \m' 



M = — shdr(fcr) = 



c 2 W 



and N = N' /c 2 . Introducing the constants a 1 = c i a 1 b 1 = 
c 4 b, e' = c 2 e, g 1 = c 2 g, ml 



(3' = cf3, and defining P' and Q' similarly to 
and Q, we get 

Then the (x, y) bloc of the metric becomes 



2 

= c z a, 
for P 



dx^_ 
~P~ 



dy 2 
Q 



, fdx' 2 



d_f_ 

Q' 



To get the appropriate transformation for the other terms 
in the metric, we have to change the coordinates as 
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t = ct' and a = c 3 a' . It follows that the two sets of con- 
stants, a, b, e, g, k, m, n, a, (3 and a' , 6', e', g', k', to', n', 
a', (3 1 describe diffeomorphic metrics. Since no other sim- 
ilar freedom exists, there are eight physical parameters 
in the charged Wahlquist metric. One can use the above 
freedom to set one of the constants to some prescribed 
value, for example one may set k = 1 or (3 = 1. The con- 
stants a and b can only be scaled by a positive constant. 
Hence they may be put equal to ±1 or 0. 
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